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稻 纵 卷 叶 蜡 种 群 系统 动态 模拟 
KAT REE 


《湖南 农学 院 > 长 沙 410128) 


摘要 本文 组 奸 了 称 纵 卷 叶 况 种 群 系统 模 似 模型 。 它 由 系统 亚 模 型 和 总 模型 二 个 部 分 组 成 。 系 统 亚 模 
型 包括 发 疹 、 死 亡 和 繁殖 三 个 子 模 块 。 系 统 总 模型 是 作 竹 提出 的 新 模型 * 它 综 含 了 前 人 所 提出 的 种 故 模 型 的 
优点 ,以 差分 方程 形式 给 出 ,以 生理 时 间 为 单位 ,考虑 了 种 群 内 个 体 间 发 育 束 率 的 差异 ,不 仅 能 模拟 种 群 数 是 
动态 ,而 且 能 模拟 种 群 年 龄 结构 ?同时 能 预测 发 和 期 。 横 型 有 效 性 检验 表明 。 模拟 结果 基本 上 能 吻合 实测 结 
果 。 文 中 还 对 5 个 影响 因子 进行 了 灵敏 性 分 析 。 


关键 词 WAER MERZA 模拟 模型 


AHIR Cnaphalocrocis medinalis Guenée 是 我 国 水 稻 主 要 害虫 之 一 ， 在 湖南 
省 每 年 发 生 4 一 5 代 ， 第 二 .四 代 分 别 为 害 早 稻 和 上 晚稻。 前 人 对 稻 纵 着 叶 蜡 进 行 过 大 最 研 
究 , 积 累 了 许多 个 体 生态 和 种 群生 态 的 资料 。 本 文 是 在 前 人 工作 基础 上 ,按照 大 系 绕 的 分 
和 解 与 综合 原则 ,探讨 种 群 系统 的 建 模 方 法 ,逐日 模拟 稻 纵 着 叶 蜡 的 种 群 数量 及 年 龄 结构 动 
态 ,为 系统 管理 与 决策 打下 基础 。 

x 据 来 源 

(一 ) 稳 纵 卷 叶 蜗 各 虫 态 有 效 积 温 及 发 育 起 点 温度 , 引 自 陈 永年 (1986)。 但 因 其 幼虫 
阶段 没有 分 龄 ,我 们 假定 幼虫 各 龄 的 发 育 起 点 温度 相同 ,幼虫 各 龄 有 效 积 温 ， 系 参照 华南 
农学 院 (1981) 关 于 第 二 代 各 龄 龄 期 资料 , 由 幼虫 期 总 积温 折算 而 来 。 

(=) 不 同 恒温 下 稳 纵 卷 时 蝇 备 里 态 存活 率 , 引 有 自 陈 永年 (1986, 内 部 资料 )。 

(=) 食料 条 件 对 黎 纵 卷 叶 嵌 各 上 虫 态 存活 率 及 成 虫 生殖 力 的 影响 , 引 自 但 建国 (1988， 
内 部 资料 )。 i 

《四 ) 长 沙 地 区 第 二 代 稻 纵 卷 叶 旦 生命 表 资 料 , 引 自 潘 桐 等 (1987) 

(CE) 温度 影响 成 虫 产 卵 三 的 资料 , 引 自 庞 雁 飞 等 (1982)。 

(A) 产 卵 时 间 分 布 ， 系 作者 观察 的 结果 ， 日 产 卵 最 占 总 产 卵 蔓 百 分 比 : 第 1 天 
40.1% ,第 2 天 34.8%, 第 3 天 15.3% ,第 4 天 9.8%o 

CE) 模型 有 效 性 检验 的 资料 ， 系 作者 1987 年 大 田 验 证 试验 的 结果 。 大 田 验 证 试验 
分 两 类 ; 一 类 是 第 四 试验 , 6 月 10 日 接 新 产 的 卵 200 粒 / 百 丛 , 接 虫 前 10 RIP HRA H+ 
施 药 清 除 其 他 病 虫 , 侯 四 面积 4X1.5 米 ?, 小 区 施 纯 氢 12.5 公斤 / 亩 ; 每 隔 5 天 调查 一 次 忠 
A LM Eiko 另 一 类 是 大 田 普 查 。 

〈《 八 》 气象 资料 来 自 湖南 省 气象 观 象 合 ,该 合 距 试验 名 2.5 公里 。 


CF 1988 年 12 月 临 到 。 
*# 国家 自然 科学 荡 金 项 目 。 研 究 过 程 中 ,得 到 本 院 陈 永年 教授 和 但 建国 同志 的 大 力 支 持 深 形 谢 误 。 
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系统 亚 模 型 
种 群 数 量 与 年 龄 结构 的 变动 由 发 育 、 存 活 ` 签 殖 , 迁 人 与 迁 出 等 五 个 子 过 程 所 控制 。 我 
们 以 稳 纵 卷 叶 虹 第 二 代 作 为 模拟 对 象 ,在 潮 南 第 二 代 是 本 地 发 生 代 , 不 存在 迁 人 与 迁 出 问 
题 ,所 以 路 去 此 琴 个 子 校 块 。 
《一 ) REFRE 
节肢 动物 种 群 有 虫 态 和 龄 期 之 分 ,在 进行 系统 模拟 时 ,我 们 定义 每 个 虫 态 为 一 个 发 育 
阶段 ,每 个 发 育 阶 段 又 等 分 为 若干 个 年 龄 级 ,以 生理 时 间 囊 示 各 虫 态 的 发 谊 历 期 ， 应 用 有 
效 积 逮 模 型 模拟 种 群 的 发 育 过 程 。 当 日 育 将 积温 的 计算 公式 如 下 : 
TCK) — COQ) TCK) 2 coc) 
BO U,K) = to T(K) < COU) C1) 
式 中 BOU.K) SBI TRAE RAR). BK AWA, TR) 为 第 下 
日 的 日 平均 温度 ，COCT)》 OS I NAAR AME 
对 于 种 群发 育 时 间 分 布 的 模拟 ,我 们 仿照 髓 车 熏 型 中 的 蕉 进 比例 法 ,假定 种 群 从 上 一 
me, 人 人 下 一 年 龄 级 时 ,六 不 是 100% 地 移动 ， 而 是 按照 一 定 的 比例 向 前 转移 ,定义 这 
个 是 俩 为 推进 比例 FOS AURIS BE EE FE TREE LE OD, Goudriaan 
(1973) 2A y HESE EENET ie At: 





sD 让 l 
p00)| C) 
式 中 FOU) AOU SUJ DOC) SAA I MEREEN ERRA, RBAN 
差 和 平均 发 育 历 期 。 发 育 历 期 方差 由 陈 永年 (1986) 的 资料 , 按 下 式 进行 计算 : 


n— i 
式 中 天 为 各 批 试 虫 实测 积温 值 ，K 为 计算 积温 值 , f OR, ?为 个 体 总 数 。 (3) 
(=) 死亡 子 模块 
在 稻 纵 卷 时 晤 种 群 系统 的 异 氟 中 ,死亡 子 横 坎 包 揪 三 个 子 槛 型 ， 即 温度 与 死亡 、 食料 
条 件 与 死亡 .天敌 与 死亡 子 模 型 。 
1. 恒温 与 死亡 ”根据 陈 永年 (1986) 资 料 ， Le RE SERIE ec EL 
成 二 次 曲线 : 


FOCI) = 1 — AOU)» | 





SOCI) = 一 126.08 + 19.047 一 0.407T? 
$O(2) = —155.47 + 20.087 — 0.409T? 
SO(3) = —208.09 + 23.97T 一 0.465T? 
SO(4) 一 一 212.18 十 24.427 — 0.4727? 
SO(5\ = —222.69 -+ 25.23T — 0.483T? (4) 
SOCb) 一 —209.62 + 24.4 一 0.468T? 
SOCI) 一 一 215.46 + 24.557 一 0.47277 
式 中 SOCLI—SO(7) 分 别 为 卵 , 1 一 5 龄 幼虫 及 师 的 存活 率 , 7 为 温度 "Do 
2. 食料 条 件 与 死亡 ”用 叶片 售 氮 量 表 示 食 料 条 件 , 报 据 但 建国 (1988) 资 料 , 食 料 条 件 


“3 其 PRPS: LALA WORD TERA Sesh SUID 313 
3 期 Yo 


对 卵 、 肾 的 存活 率 影响 不 显著 ,幼虫 的 存活 率 则 与 叶片 合 氮 量 呈 直线 关系 : 
SN(2) = 37.26 + 10.94N (r = 0.98) 6) 
SN(5) = 65.47-+7.42N (r= 0.93) 

式 中 SN(2),SN(5) 分 别 为 低龄 (1 一 3 AEA C45 I BRRR , NE Se 

Ait. | 

3. 天敌 作 用 后 存活 率 ， 天 敌对 害虫 的 捕食 或 寄生 作用 ,关系 复杂 ,目前 尚 缺乏 实用 性 
较 强 的 数学 模型 ,因此 我 们 引用 潘 桐 等 (1987 ) 的 生命 表 资 料 , 将 各 各 天敌 因子 作用 后 的 存 
活 率 累 乘 ,得 到 天 敌 作用 后 的 总 存活 率 。 

4. 自然 条 件 下 各 年 龄 级 存活 率 的 模拟 。 上述 子 模型 都 是 表达 各 种 因子 对 阶段 存活 
率 的 影响 。 为 了 便于 系统 模拟 ， 还 需 建立 各 年 龄 级 存活 率 子 模型 。 假 定 同 一 发 育 阶 段 各 
年 龄 级 具有 相同 存活 率 ( 陈 维 博 等 ,1981), 则 有 : 

SRU) = 一 -一 一 一 一 (6) 
FOU 
1 FOC) + 400), RCI 
AH SRC), RU) 分 别 为 年 龄 级 和 发 育 阶 段 的 存活 率 ，FOCI) 为 第 了 阶段 推进 比例 ， 
4071) 为 了 阶段 年 龄 级 数 。 

为 了 模拟 自然 条 件 下 各 阶段 总 的 存活 率 ， 假定 温度 和 食料 条 件 对 存活 率 的 影响 不 存 
在 时 汪 作 用 ,因而 第 了 工 阶段 、 第 到 日 条 件 下 的 阶段 总 存活 率 〈8ST(7, 有 7)) 为 : 

ST(1,K) = SOCI) » SN'(I) < SP(1) X 107 (7) 
式 中 SPC) 为 天 敌 作用 后 的 存活 率 ; SOU) 为 温度 影响 的 存活 率 ;出 (4) 式 求 出 ， 计 算 
尘 (4) 式 中 的 TT 取 第 及 日 平均 温度 ; SNU) 为 食料 条 件 影 响 的 存活 率 , 由 (5) 式 将 低龄 幼 

虫 和 高 龄 幼虫 的 存活 率 等 分 给 幼虫 各 岭 B: 
SN' (2) = SN'(3) = SN'(4) = W SNC) (8) 
SN'(5) = SN'(6) = W SNCS) 

SN'(2)—SN'(6) 分 别 为 1 一 5 龄 幼虫 食料 条 件 影 响 的 存活 率 。 
(=) 策 殖 过 程 子 模块 
1. 温度 与 繁殖 力 ” 适 家 温度 下 恒温 与 成 虫 产 卵 最 的 关系 模型 引 自 庞 雄 飞 等 (1982); 
412.45 一 12.06T T>26%° - (9) 


100 
1 + exp(10.07—0.56T) 


式 中 RT AFKEER, TAEC 假定 变温 环境 对 产 卵 生 的 影响 不 存在 时 滞 
效应 , 则 可 将 (9) 式 中 的 温度 了 代 之 以 日 平均 温度 T(KE)， 从 而 模拟 玉 日 变温 条 件 下 成 虫 
产 卵 熏 。 
2. 食料 条 件 与 繁殖 力 ”食料 条 件 以 施肥 水 平 表示 , 并 以 每 亩 12.5 ATARE RR 
准 施 肥水 平 ,定义 孵 授 达标 准 施肥 水 平 的 百分率 RN ,为 任 一 施肥 水 平 下 产 卵 下 对 标准 
施肥 水 平 下 产 卵 重 的 百分比 。 根 据 但 建国 (1988) 的 资料 可 拟 合成 曲线 方程 : 
RN = 54.12 十 3.04F, 十 5.4277F3 (10) 





er ={ 
T & 26% 
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其 中 FP, 为 施肥 水 平 。 
3. 逐日 产 卵 量 的 模 氟 ”首先 根据 作者 观察 资料 ， 将 第 天 日 产 卵 景 占 总 产 卵 是 的 百 分 
比 GCK), 对 产 卵 日 期 (开始 产 卵 后 的 天 数 ) 氢 合成 产 旷 时 间 分 布 模 型 : 


G(K) = 14.44 + 22.01K — 6.1K? (11> 
BAO) RE RSA PRE PON (RE): 
R(K) = RT + RN + G(K)/100 (12) 


系统 总 模型 


节肢 动物 自然 种 群 不 仅 存在 年 龄 结构 ， 和 而 卫 同 年 龄 组 个 体 间 的 发 育 速 率 也 不 一 致 
LesLie (1945) 提出 了 包含 年 龄 结构 信息 的 算 阵 位 型 , RP HEAR (1981) 建立 了 适用 于 变 
温 杂 件 下 模拟 的 变 维 矩阵 异型 ，Van Foester (1959) 提出 了 包含 年 龄 结构 的 偏 筱 分 方程 
模型 , 但 上 述 语 焉 吕 没 有 著 感 个 体 阅 发 育 历 期 的 符 异 。Goudriaan (1973) FEW AS E 
MEERI Cbex-car model), MIRRE HNP RAEI o Ruesink (1976) Ht 
BEY 6 EI SUS US ME op Dy A ee 2 A r ET r Pb Ok ri HEE 
HUE DZ PEL BITRE ESE WA ee Bb BD 3 BE Ey ME Be A ed Ae 
HEIER A LA He de ih RE HEU Pe PT BET yg PE ile 3 
FES eS TEAC SRA BILA THE HS Se ER DIRT o Mie PHL SB fe ae 
BLS 28 oF PROP PDS RT E Fo E TE EE DR STE 
FRR 

(一 ) 定义 与 说 明 

1. ASE PI RSS Bie ike RID RNP RAE Bor LEAP 
RARE BTM RMEERINE 5) 400) 个 年 龄 等 级 。 

DOU 
AOC) = poU) (13) 





式 中 DOC) ABI MEAAG AT 为 年 龄 等 级 的 期 中, 并 以 AT 作为 模拟 步 长 。 
2. 第 7 Mi BEM AES K KITEE EEOAE A 
m BO(I,K) 
NS (I,K) INT ( Ue +05 ) (14) 


式 中 BOU,K) 35 1 REBR HARNI, HORRET R, INT ARERR 
3. 梳 拟 时 间 每 移动 一 个 步 长 ， 了 阶段 因 种 岩 从 上 一 年 龄 狼 发 育 到 下 一 年 龄 级 向 扒 进 
tha FOC), FERH 8RCDD， 而 停留 于 原 年 龄 级 不 动 的 比例 为 1 一 FOU). #8 
率 为 100%, 
4. 由 于 各 阶 钻 发 育 起 点 温度 不 相同 , 故 同一 天 不 同 发 育 阶段 只 育 不 同 的 有 将 积温 值 。 
当 在 第 玉 天 内 发 生 阶 段 之 间 的 转移 有 时, 如果 已 名 种 群 订 第 了 阶段 内 猕 移动 了 个 年 龄 级 , 那 
人 么 当日 肉 在 第 1-1 阶段 内 可 移动 的 最 大 年 龄 级 数 VU 一 1) 为 : 
VU 一 中 一 [N8CR) 一 中 .BOC —1,K)/BOCI,K) (1s) 
(=) 异型 的 差分 方程 表示 式 
X(K,1,1) = [1 — FOC) IM . X(K — 1,1,1) 
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AQ(N) 
+ 5} 0.5+RCK, M)X(K —1,N,M) (16) 


Mm 1 


J 
X(K,1,D= >> [BROCDJ7T 


M= J NB K) : 
© [1 — FOCI) MOH. [SROV 
. XK = 1,1,M) . 

J>N8U,K) +1 (17) 


J 
XK, 1,0) = STL FOC) P” + [1 — FOC) ns 。 
M=1 


[SRCI + XCK— 1,1,M) 
AOA) 


十 > [FOCI — 1)]400-9-M 


Mn dQUI—1)—-Vd—1) 
e [FODDI . [1 一 FOCI -一 1) ]¥-400-9+ van) 
e X(K— 1,1 —1,M) + [SRCI [SRC — 1)] 40079 
J<N8(I,K) ~ (18) 

式 中 X(K,1,J) EREKE, BINE, BI 年龄 级 的 数量 ，RCK,M) HRKABM 
年 龄 级 成 虫 的 产 卵 量 , 其 它 符号 含义 见 前 述 。(16) 式 计算 第 天 日 卵 期 第 一 年 龄 级 的 数量 。 
(17) 式 计算 第 玉 日 .第 工 阶段 .大 于 当日 最 大 推进 年 龄 级 数 的 那些 年 龄 级 的 数量 ， 它 们 人 金 
部 由 本 阶段 内 低 年 龄 级 发 育 及 本 年 龄 级 滞留 而 来 。(18) 式 计算 第 KK 日 ,第 了 阶段 小 于 或 
等 于 当日 最 大 推进 年 龄 级 数 的 那些 年 龄 级 的 数量 ， 或 者 由 本 阶段 低 年 龄 级 发 育 及 本 年 龄 
级 滞留 而 来 ,或 者 由 前 一 阶 眉 发 育 而 来 , 意 即 不 仅 发 生 阶 跟 内 的 转移 ， 同 时 发 生 阶段 间 的 
转移 。 


计算 机 模拟 及 有 效 性 检验 | 
在 计算 机 模拟 时 ,将 稻 纵 卷 叶 旺 的 一 个 世代 分 成 "个 发 育 阶段 ( 卵 、1 一 5 hh 
产 卵 前 期 、 产 卵 期 )， 年 龄 级 期 距 为 8 个 日 度 。 根 据 陈 永年 (1986》、 潘 桐 等 (1987) 资 料 , 得 
出 各 阶段 年 龄 级 数 及 天 敌 作用 后 存活 率 ， 列 于 表 1 。 模 氟 时 ， 计 算 机 内 存 人 表 中 全 部 次 
Ph 运行 时 输入 移 栽 后 天 数 、 施 肥水 平 、 日 平均 温度 ,种群 的 初始 状态 (年 龄 结构 及 数量 ); 


Rl MAST RAMRERRASER A as 



























阶段 序号 8 9 
相应 的 下 六 NTT. 产 琵 其 
有 效 积温 (日 度 ) 32.19 | 32.19 | 53.65 | 86.30 | 21.72 | 22.00 
STIG AA BECO) | 16.23 | 16.23, ‘| 19.00 | 19.00 
年 龄 级 数 3 | 3 
推进 比例 0.938 0.916 | 0.916 | 0.916 | 0.916 | 0.916 | 0.936 | 0.923 | 1.000 
~ EE a a 
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结束 时 输出 逐日 年 龄 结构 及 种 群 数量 。 模 拟 程序 用 BASICA 写成 , 在 IBM PC/XT 机 
上 实现 。 

为 了 检验 模型 的 有 效 性 , 我 们 输入 施肥 水 平 12.5 公斤 / 亩 纯 氮 , 移 裁 后 38 KR 
是 为 200 粒 / 百 从 卵 。 模 拟 结果 与 实测 结果 基本 上 接近 ， 说 明 模 型 有 效 性 较 好 (了 见 表 2)。 
但 1 龄 幼虫 实测 结果 严重 偏 低 ， 主 要 原因 可 能 是 1 龄 幼虫 难以 调查 准确 。 表 中 所 有 实测 
值 痢 有 偏 低 的 趋势 ， 主 要 因为 实测 值 是 笼罩 试验 的 结果 ， 第 四 条 件 与 一 般 大 田 条 件 有 差 
异 。 另 外 大 田立 查 结果 也 与 模 掀 结果 基本 上 吻合 ,说 明 模型 具有 一 定 的 实用 价值 。 
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S oO 5 fa] S o 3 oO RY (2) 
38 200 200 
42 175.11 | 165.0 
47 45.8 23.0 92.17 ata 
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He: “RRRA ats 5.0 3} SRR US ilo 


RM DA 

PRR ee ES, API REER ARER KE Se F 
BINS SY BIBER — 20 % 10% H107 20 %o FEL F HOPE ZK AF DS CE he ME) 
He SF EE Ss FAS A, MA TRA RT TAF SS (eM AY TEBE, AURAL 3 。 
由 点 3 TA ERREN SRT BE REER i SS HB 9 E 
SS WOMB ASA RZ 3 AOA ERASE Me FUSE APRS AS APEC US, ph ee AE BSE 
BEN RBIS ESB BORO HA Se BEA, SRBC E HERR AG REEKSE Be Be 5] 
HARP APES TABU aR BE SWUM PO PNUD BR I 38 R, IRE AREA 
时 ,种 铬 增长 最 快 。 日 平均 温度 的 升 高 , 比 日 平均 温度 的 隆 陶 ， 对 种 群 的 抑制 作用 更 大 。 
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害虫 发 生 期 









日 平均 温度 





-20 一 29.0 —30.1 一 56.2 0 
一 10 一 17.0 一 16.5 一 45.3 0 
+10 一 36.1 一 56.4 十 34.2 0 
+20 0 






讨 论 


《一 ) 本 文 提出 了 模拟 节肢 动物 种 群 动态 的 一 个 新 模型 ， 它 综合 了 前 人 所 提出 的 种 群 
动态 模型 的 优点 ， 以 生理 时 间 为 单位 ， 能 够 横 拟 个 体 间 发 谊 速率 的 差异 和 种 群 的 年 龄 结 
构 。 它 以 差分 方程 的 形式 给 出 ,这样 便 子 在 计算 机 上 运行 。 模 型 有 效 性 检验 表明 ,模拟 结 
果 基 本 上 能 吻合 实测 结果 ,这 说 明了 模型 的 合理 性 和 可 行 性 。 

(=) 文中 组 建 稻 纵 卷 叶 蝇 种 群 系统 亚 模型 时 , 对 于 自然 条 件 下 种 群 存活 率 和 繁殖 力 
的 模 氢 ,是 假定 温度 和 食料 条 件 等 因子 对 存活 率 和 繁殖 力 的 影响 不 存在 时 涪 作 用 ,通过 恒 
定 条 件 下 存活 率 和 繁殖 力 模拟 第 玉 天 条 件 下 的 存活 率 和 繁殖 力 ， 这 里 隐 含 着 把 1 天 中 的 
温度 和 营养 条 件 视 为 恒定 的 假设 ,这 种 做 法 具有 一 定 的 代表 性 。 但 时 澡 作用 的 考虑 ;有待 
进一步 研究 。 

(=) Bae TR A BR, 有 得 深入 研究 ,以 便 完善 种 群 
系统 模型 。 
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SIMULATION OF THE DYNAMICS OF POPULATION SYSTEM OF THE 
RICE LEAF ROLLER, CNAPHALOCROCIS MEDINALIS GUENEE 


ZHANG JUN-PING CHEN CHANG-MIN 


(Hunan Agricultural College, Changsha 410128) 


The general model and its submodels of population system of the rice leaf roller Cnapha- 
locrocis medinalis are developed according tothe principle of decomposition and synthesis of 
huge system. The system submodels consist of three parts: developmental submodel, mortality 
submodel and reproduction submodel. The general model of population system is a new kind 
of difference equation model which is suggested first in this paper. Our new model synthes- 
izes the advantages of population simulation models which have been suggested ever before. It 
uses physiological time as the unit of time to simulate the developmental process. The var- 
dance of developmental rate among individuals has been considered th the model. It can be 
used to simulate not only experimental population, but also natural population. It can simu- 
late the dynamics of population density, age stracture and developmental stage day by day at 
the same time. The simulation results can approach the observed results satisfactorily. he 
sensitivity ahalysis on five factors is conducted in the paper. 
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